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ABSTRACT 

We have carried-out 98-level configuration-interaction / close-coupling (CI/CC) inter¬ 
mediate coupling frame transformation (ICFT) and Breit-Pauli /^-matrix calculations for the 
electron-impact excitation of Be-like A1 The close agreement that we find between the 
two sets of effective collision strengths demonstrates the continued robustness of the ICFT 
method. On the other hand, a comparison of this data with previous 238-level CI/CC ICFT 
effective collision strengths shows that the results for excitation up to « = 4 levels are sys¬ 
tematically and increasingly underestimated over a wide range of temperatures by /^-matrix 
calculations whose close-coupling expansion extends only to « = 4 (98-levels). 

Thus, we find to be false a recent conjecture that the ICFT approach may not be com¬ 
pletely robust. The conjecture was based upon a comparison of 98-level CI/CC Dirac /^-matrix 
effective collision strengths for A1 with those from the 238-level CI/CC ICFT /^-matrix 
calculations. The disagreement found recently is due to a lack of convergence of the close¬ 
coupling expansion in the 98-level CI/CC Dirac work. The earlier 238-level CI/CC ICFT work 
has a superior target to the 98-level CI/CC Dirac one and provides more accurate atomic data. 

Similar considerations need to be made for other Be-like ions and for other sequences. 

Key words: Atomic data - Techniques; spectroscopic 


1 INTRODUCTION 


Electron-impact excitation is the dominant process for populating 
the radiating states of ions whose emission lines form the basis for 
the spectroscopic diagnostic modelling of non-equilibrium astro- 
physical and laboratory plasmas. As such, a great deal of effort 
over many years has gone in to calculating a large amount of colli¬ 
sion data, which is incorporated into databases and modelling suites 
such as CHIANtB and OPEN-ADA;^. The pre-eminent method¬ 
ology used is the f?-matrix one. However, there are many variations 
on this theme. 

The intermediate coupling frame transformation (ICFT) R- 
matrix method is an approximation to the Breit-Pauli /^-matrix 
method (BPRM) which neglects the spin-orbit interaction between 
the colliding electron and the ion. Thus, it is physically well mo¬ 
tivated. There is a good deal of literature which verifies the ac¬ 
curacy of the ICFT approach: the original compariso ns of the re¬ 
sults of ICFT and Breit-Pauli R- matrix calculations by Griffin et"^ 
( Il998lll999h on Mg- like ions andiBadnell & GriffinI l l 19991) on NF^: 
more recent ones bv lUiang & Badn^ 1 I 2 OIOI) comparing ICFT R- 
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matrix and Dirac Atomic R-matrix Code (P ARC) calculations for 
Ne-like Fe*'^^ and Kr-'^+ ■ lUiang e7 iD ( l2009h on ICFT with Breit- 
Pauli andDARC_forNa4ike^e^^^j_^nd_mos^recentl^Badnell & 
Ballance l l20l4) compared the results of ICFT. BPRM and DARC 
calculations for Fe^^. The differences observed between ICFT and 
other R-matrix results are all well within the uncertainties to be 
expected due to the use of different configuration interaction and 
close-coupling expansions and resonance resolution. Indeed, the 
iBadnell & Ballanc3 j2014l) work on Fe^^ used identical atomic 
structures and close-coupling expansions for the ICFT and Breit- 
Pauli R-matrix calculations and found excellent agreement, better 
than 5%. Their DARC calculations out of necessity used a some¬ 
what different atomic structure, and the low-level structure of Fe^^ 
is challenging, but still gave agreement to within ~ 20% at 10“^ K. 

Extensive calculations have been carried-out applying the 
ICFT R-matrix method to whole isoelectronic sequences (for el- 
ements, typically, up to Zn), the m ost recent ones being Mg-like 
jperntodez-Menchero et alj2014bl). Be4ike_(Fernandez-Menchero 
et al. l2014al) and B-like jUiang et alj|2012h — this last paper con¬ 
tains references to earlier sequences. 

Thus it is both surpris ing and of great conce r n to fin d a work 
which counters this trend: lAggarwal & KeenanI FoiSah made a 
comparison of the results of a 98-level D ARC calculation on Be¬ 
like A1 l lAggarwal & Keenarj|2014ch with a 2 38-level ICFT 
R-matrixone_jFemandez;Mencheio_et_ffi 1M)- _Aggarwal & 
Keenan T20I5ah found that the effective collision strengths obtained 
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from the 238-level ICFT f?-matrix calculation were significantly 
larger than the P ARC ones in many instances and th ey suggested 
that the results of iFernandez-Menchero et ^ ( 12014all were less re¬ 
liable, querying the ICFT method, resonance resolution and its high 
energy / tempera ture behaviour. _ 

In addition, lAggarwal & Keen^ (l2015al) alighted on the re¬ 
cent paper by IStorev et alj ( |2014|) which reported a problem in 
the outer-region ICFT calculation of when compared to a full 
Breit-Pauli calculation, and suggested that thi s could be the main 
cause of the discrepancies for Al . However, IStorev et alJ ( l2014li 
noted that such an issue only arises when resonance effective quan¬ 
tum numbers become small. The problem is peculiar to low-charge 
ions such_^sO^^_andjinusuall^_small_R;niatrix__box_size^Storey 
et al. ( |2014) focused on providing a solution to their problem at 
hand. We note that the problem does not arise in the first place 
if the R-matrix box size is increased (beyond its default in their 
case) to encompass a spectroscopic n = 3 orbital, say. This is why 
the issue had not arisen before: all previous calculations (including 
iFemandez-Menchero et alj i2014ah) used larger, ofte n much larger, 
box sizes. The problem noted by Store^rt is not relevant 

in general. lAggarwal & KeenanI ( l2015all noted a similar trend for 
other ions_oftheBedike_se 2 uenceCl^^^jJC^^_|jGe^^^Aggarwal & 
Keenan l2014ali and Ti‘^^ dAggarwal & Keenanll2012h . 

Whe re does that leave us with re gard to the discrepancies 
noted by lAggarwal & KeenanI (l2015al)? The concern of Aggar- 
wal & Keenan j2015ah was the disagreement between R-matrix 
calculati ons of (apparent) comparable complexity. However, the 
worksofjA^garwal_&_^eenaijj2014^^_an^Fernandez-Menchero 
et al. l l2014ah are not of comparable complexity. We show here 
that the much lMger_close^cgygfing_ex^nsion_usedb^Femandez- 
Menchero et al. l2014all (238 vs 98 levels) gives rise to a system¬ 
atic enhancement of effective collision strengths over a wide range 
of temperatures, which increases as one excites higher-and-higher 
levels. 

In addition, we analyze the uncertainty in the effective colli¬ 
sion strengths due to the incompleteness of the configuration inter¬ 
action expansion, the validity of the ICFT vs Breit-Pauli method, 
viz. the neglect of the spin-orbit interaction of the colliding elec¬ 
tron, and the effect of resonance resolution and position on low 
temperature effective collision strengths. 

The paper is organized as follows. In section[2]we describe the 
methodology we used for the different calculations we have per¬ 
formed. In section [3] we discuss the atomic structure of Al and 
present results for energies, line strengths and infinite energy plane 
wave Born collision strengths. In section |4] we compare and con¬ 
trast effective collision strengths. In section|5]we present our main 
conclusions. Atomic units are used unless otherwise specified. 


2 METHODOLOGY 

In the following sections we compare the results of the 238-level 
configuration-i nteraction / close-coupling (CI/CC) ICFT R-matrix 
calculation by iFernandez-Menchero et al] ( 12014al) for Al with 
the results of new ICFT and Breit-Pauli 98-level CI/CC R-matrix 
calculations — t he latter being the same sized Cl & CC expansions 
that were used bv lAggarwal & KeenanI ( l2012Ll2014allcli . Where pos¬ 
sible (meaningful), our new 98-level CC calculations f ollow the 
same prescription as IFernandez-Menchero et alJ(l2014jl e.g. with 
respect to angular momentum and energy specification. 

Thetarget^escrigtioiruses^lKAiTOSTRUCTURE^rograin^Bad- 
nell lTOllh . IFernandez-Menchero et a lJ( l2014ti) included all of the 


Table 1. Scaling paiameters (A„i) used in the 98-level Cl stmcture calcula¬ 
tion. 


Orbital 

A„i 

Is 

1.57175 

2s 

1.29675 

2p 

1.17081 

3s 

1.29004 

3p 

1.15958 

3d 

1.29791 

4s 

1.29797 

4p 

1.15413 

4d 

1.30212 

4f 

1.50916 


configurations Is^ {2s^, 2s2p, 2p-j and Is^ {2s, 2pj nl with « = 3 - 7 
and / = s, p, d, f, g for n = 3 - 5 and / = s, p, d for n = 6,7, which 
makes a total of 238 levels. In the present work, we restrict our¬ 
selves to n = 4 (with Z = s, p, d, f) which gives a total of 98 levels. 
The Thomas-Fermi potential scaling parameters, T,,;, for the atomic 
structure calculation with the Cl basis set of 98 levels are given 
in Table [J those for tbe Cl calculatio n with 238 levels are given 
in IFernandez-Menchero et^ (l2014al) . The T,,; parameters for the 
98-level Cl calculation were obtained in the same way as the one s 
for the 238-level calculation dFemandez-Menchero et aul2014ah : 
by minimizing the equally-weighted sum of all LS -coupling term 
energies. 

For the collision calcula t ions w e use the inner-region R-m atrix 
programs of iHummer et alj (Il993h: iB errington et al.l(l l995l) and 


the ou ter-region stop program of Bemngtonrt^ 1 19871) : Badnelll 


d 19991) plus the ICFT one o f iGriffin et alj dl998l) . We will com¬ 
pare the results of ICFT and Breit-Pauli R-matrix methods to treat 
relativistic effects in the scattering problem. Both can use the ex¬ 
act same Breit-Pauli atomic structure. This is important as it en¬ 
ables us to isolate differences due solely to the di ffering treatment 
of rela tivistic effects in the scattering. We note that iBerrington et aP 
d2005l) have shown that the Breit-Pauli R-matrix method gives es¬ 
sentially the same results as the Dirac one for Z < 30. 

The ICFT R-matrix method first carries-out an L5-coupling 
close-coupling calculation for the Cl target described above and it 
includes the mass-velocity and Darwin one-body relativistic oper¬ 
ators. The LS -coupling reactance R-matrices are first recoupled to 
/R-coupling and then transfor med to intermediate coupling using 
the term coupling coefficients dHummer et al.|[T99^ for the Breit- 
Pauli target. This imposes the exact same atomic structure on the 
ICFT R-matrices as a Breit-Pauli one which uses the same Cl ex¬ 
pansion and radial orbitals. In particular, the fine-structure levels 
within a term are non-degenerate in the final scattering calculation. 
This meftod^s^^e_^ne_^Wch_^as_^se^^^Fernandez-Menchero 
et al. d2014al) . 

The second method which we use is the Breit-Pauli (BP) one. 
In this method the close-coupling expansion is made in intermedi¬ 
ate coupling as well, in addition to the target Cl, i.e. the one-body 
effective (nuclear plus Blume & Watson) spin-orbit operator is in¬ 
cluded in the (N + l)-electron Hamiltonian. The Breit-Pauli for¬ 
malism increases considerably the size of the Hamiltonian matrix, 
which makes it impractical to use this method for 238 CC levels. 

The essential physical difference between the ICFT and Breit- 
Pauli R-matrix methods is the neglect by the former of the effective 
spin-orbit interaction between the colliding electron and the ion. 
The practical benefit of the ICFT method over the Breit-Pauli one 
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Validity of the ICFT R-matrix method 3 


is the diagonalization of much smaller (A1+ l)-electron Hamiltonian 
matrices and a much smaller set of coupled scattering equations to 
he solved in the outer-region hy stgf. 


3 STRUCTURE 


Table[^compares our energies for the first 98 levels of A1 calcu¬ 
lated with the 98- and 238-level Cl targets. These energies are com¬ 
pared also with the observed one s tabulated in the nistH database, 
which were taken from work of iMartin & Zalub^ l ll979ll . There 
are ten levels in Table[2]which do not follow the same order in both 
structure calculations: from level index 82 to 84 and from 88 to 
94. We use the order of the 98-level Cl calculation to index levels 
for comparison purposes. Energies calculated with both Cl expan¬ 
sions have differences smaller than 0.5% with the observed ones for 
most of the levels, and differences of around 2 % for the low lying 
si nglet states. In general , the dif ferences are similar to those found 
bv lAggarwal_&^Keenan| j2()j_4j}HnJhej)riginalj)aperfrom^ Martin 
& Zalubas ( Il979l) there are several gaps in the level energies, and 
some of them are labeled as “inaccurate”, so we refer the reader to 
the original work to avoid hasty conclusions. 

It is difficult to relate differences in energies directly to dif¬ 
ferences in collision data. However, differences in oscillator or 
line strengths (5) and infinite energy plane-wave Bom collision 
strengths (fl™®), essentia lly non-dipole e l ectric multipole line 
strengths, can be so related. lSurgess & Tullvl l l 19921) show how infi¬ 
nite energy / temperature and ordinary / effective collision strengths 
(Oco/Too) from any scattering calculation, including an /^-matrix 
one, are directly related to these quantities, viz. 

( 1 ) 


for non-dipole allowed transitions, while for electric dipole ones 




4S 

/ E 

Ooo = 

lim 

—In 



£-»oo 

3 

[ae 

and 


4S 1 

ikT 

T„„ = 

lim 

—In 



r->oo 

3 1 

[AE 


( 2 ) 

(3) 


where A£ is the excitation energy for the transition. 

In practice, we find that changes in the line strength, S, or 
n™® between two different atomic structures change not only the 
infinite energy values but also the (background) ordinary collision 
strength correspondingly over a wide range of collision energies, 
and hence the effective collision strength over a wide range of tem¬ 
peratures, unless dominated by resonances. A 20% change, say, in 
S or provides a very realistic measure of the resulting change 
in effective collision strengths. Thus, care must be taken when at¬ 
tempting to deduce anything about scattering methods from differ¬ 
ences in the collision data without reference to differences in fhe 
underlying atomic structure. 

To illustrate, we show in Fig.[T]a comparison of the reduced 
quantity ;/„, given by 

= or (4) 

for the two atomic structure calculations which we consider (which 
we label 98-level Cl and 238-level Cl). We show points for a total of 
4035 transitions resulting from the common 98 levels. These split 


^ http://physics.nist.gov 



Figure 1. A comparison of yoo (see eq.|4) for the 98- vs 238-level Cl atomic 
structures for transitions amongst the 98 lowest common levels of A1 
o: transitions with upper level with n = 2; □: transitions with upper level 
with n = 3; o: transitions with upper level with n = 4; dashed lines: 20% 
fractional difference. 


into 1466 dipole transitions and 2569 Bom-allowed transitions, 
while there are an additional 718 forbidden transitions that are not 
represented. We highlight by colour and symbol transitions to upper 
levels with n = 2, 3,4. We see that transitions up to n = 2 are very 
well converged. Only one transition differs by more fhan 10%; it is 
the very weak (~ 10“**) Born transition 6-9: 2p^ ^Pq - 2p^ *D 2 (off 
the scale). There is good convergence for transitions up to u = 3, 
with 233 out of 801 transitions which differ by more fhan 20% buf 
mosfly for fhe weaker transitions. However, there is much more 
spread for transitions up to n = 4, 1834 of 3202 transitions differ 
by more than 20%. This plot illustrates that the n = A levels of the 
98-level Cl structure are not so well converged with respect to the 
Cl expansion. We expect that their representation by the 238-level 
Cl expansion to be much better converged. 

We circle several transitions for comment. The transition 31- 
41: 2p3p^Pi - 2 p 3 d^P 2 is a dipole one, its line strength changes 
by a factor of 18 between the 98- and 238-level structures. Note 
also the transitions 75 - 80: 2p4d *D“ - 2p4f'F 3 and 73 - 81: 
2p4d *’F 2 - 2p4f ^F 3 . Both transitions are dipole allowed and quite 
strong and yet differ by about a factor of 20 and 30 between the two 
structures. These levels lie towards the upper end of the 98-level Cl 
expansion. These differences highlight the need for as accurate an 
atomic description as possible to obtain the best target for accurate 
collision data and the need to exercise extreme caution when mak¬ 
ing comparisons of collision data based on different atomic struc¬ 
tures. 


4 COLLISIONS 

In this section, we carry-out a series of comparisons of effective 
collision strengths for A1 at the t emperature of peak abundance 
for an electron collisional plasma l lBrvans et al.ll20()^ . We then 
look at issues relating in particular to low temperature (e.g. pho- 
toionized) plasmas and much higher temperature plasmas (e.g. so¬ 
lar flares). 
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Table 2. target levels. Key: t level index; Conf.: configuration; Level: level IC designation (largest weight); -Enist^ observed energy from the NIST 
database iMartin & Zalubaill979l) ; E 98 : calculated energy with 98-level Cl expansion; £ 238 ^ calculated energy with 238-level Cl expansion (see text); %: 
percentage difference between theoretical and NIST data. All energies are in cm“’. 


i 

Conf. Level 

Enist 

E 98 ( 

%) 

E 238 ( %) 

i 

Conf. Level 

Enist 

Egg ( %) 

E 238 ( 

%) 

1 

2s2 

‘So 

0. 

o.( 

0) 

0. ( 0) 

50 

2s 4p 

3po 

- 

2504724. ( 

-) 

2500556. ( 

-) 

2 

2s 2p 

3po 

0 

155148. 

155722. ( 

0.4) 

155539. ( 0.3) 

51 

2s 4p 

3po 


2505125.( 

-) 

2500975.( 

-) 

3 

2s 2p 

3po 

156798. 

157487.( 

0.4) 

157404. ( 0.4) 

52 

2s 4p 

Ipo 

- 

2508073. ( 

-) 

2503514. ( 

-) 

4 

2s 2p 

3po 

160429. 

161146. ( 

0.4) 

161278. ( 0.5) 

53 

2s 4d 

^Di 

- 

2520502. ( 

-) 

2516248. ( 

-) 

5 

2s 2p 

Ipo 

^1 

300490. 

309273. ( 

2.9) 

307209. ( 2.2) 

54 

2s 4d 

3D2 


2520573. ( 

-) 

2516321.( 

-) 

6 

2p^ 

'"Po 

404574. 

408026. ( 

0.9) 

407826. ( 0.8) 

55 

2s 4d 

3D3 

- 

2520682. ( 

-) 

2516434. ( 

-) 

7 

2p^ 

'"Pi 

406517. 

409969.( 

0.8) 

409888. ( 0.8) 

56 

2s 4d 

‘D2 

2527560. 

2530218. (0.1) 

2525880. (- 

-0.1) 

8 

2p- 

^P2 

409690. 

413420. ( 

0.9) 

413526. ( 0.9) 

57 

2s 4f 

3po 

2528570. 

2530774. (0.1) 

2526258. (- 

-0.1) 

9 

2p- 

‘D 2 

449732. 

458157.( 

1.9) 

457831. ( 1.8) 

58 

2s 4f 

3po 

^3 

2528570. 

2530815. (0.1) 

2526299. (- 

-0.1) 

10 

2p- 

‘So 

553783. 

567794. ( 

2.5) 

567267. ( 2.4) 

59 

2s 4f 

3po 

2528570. 

2530871. (0.1) 

2526355. (- 

-0.1) 

11 

2s 3s 

-'Si 

1855760. 

1856089. ( 

0.0) 

1852844. (-0.2) 

60 

2s 4f 

Ipo 

- 

2533730. ( 

-) 

2528991. ( 

-) 

12 

2s 3s 

‘So 

1884420. 

1886214.( 

0.1) 

1882216. (-0.1) 

61 

2p4s 

3po 

^0 

- 

2661563. ( 

-) 

2655681.( 

-) 

13 

2s 3p 

Ipo 

^1 

1923850. 

1925826. ( 

0.1) 

1922358. (-0.1) 

62 

2p4s 

3po 

- 

2662724.( 

-) 

2656820.( 

-) 

14 

2s 3p 

3po 

^0 

- 

1928630. ( 

-) 

1925009.( -) 

63 

2p4s 

3po 

^2 

- 

2666840.( 

-) 

2661460.( 

-) 

15 

2s 3p 

3po 

- 

1929220.( 

-) 

1925611.( -) 

64 

2p4s 

Ipo 

^1 

- 

2673641.( 

-) 

2666939.( 

-) 

16 

2s 3p 

3po 

- 

1930059.( 

-) 

1926462. ( -) 

65 

2p4p 

‘Pi 

- 

2681281. ( 

-) 

2675347.( 

-) 

17 

2s 3d 

^D, 

1965860. 

1967770.( 

0.1) 

1964163. (-0.1) 

66 

2p4p 

-'Di 

- 

2684513.( 

-) 

2678722. ( 

-) 

18 

2s 3d 

3D2 

1966080. 

1967980.( 

0.1) 

1964378. (-0.1) 

67 

2p4p 

3D2 

- 

2684938.( 

-) 

2679129.( 

-) 

19 

2s 3d 

^D3 

1966300. 

1968296.( 

0.1) 

1964701. (-0.1) 

68 

2p4p 

3D3 

- 

2688348.( 

-) 

2682876.( 

-) 

20 

2s 3d 

‘D2 

1992340. 

1997586.( 

0.3) 

1993399. ( 0.1) 

69 

2p4p 

^Si 

- 

2691166. ( 

-) 

2684805. ( 

-) 

21 

2p3s 

3po 

^0 

2057140. 

2055664. (- 

-0.1) 

2050249. (-0.3) 

70 

2p4p 

^Po 

- 

2691420.( 

-) 

2685788.( 

-) 

22 

2p3s 

3po 

2057140. 

2057249. ( 

0.0) 

2051958. (-0.3) 

71 

2p4p 

^Pi 

- 

2694295.( 

-) 

2688609.( 

-) 

23 

2p3s 

3po 

^2 

2057140. 

2060910.( 

0.2) 

2055983. (-0.1) 

72 

2p4p 

^P2 

- 

2694733.( 

-) 

2689496.( 

-) 

24 

2p3s 

Ipo 

^1 

2091870. 

2090063. (- 

-0.1) 

2084057. (-0.4) 

73 

2p 4d 

3po 

^2 

- 

2697404. ( 

-) 

2691501. ( 

-) 

25 

2p3p 

‘Pi 

2094820. 

2097317. ( 

0.1) 

2093229. (-0.1) 

74 

2p 4d 

3po 

^3 

- 

2699888.( 

-) 

2694096.( 

-) 

26 

2p3p 

^D, 

2102330. 

2105510. ( 

0.2) 

2101249. (-0.1) 

75 

2p 4d 

Ido 

- 

2700756.( 

-) 

2695172. ( 

-) 

27 

2p3p 

3 D 2 

2103900. 

2107151. ( 

0.2) 

2102917. (-0.0) 

76 

2p4p 

‘D 2 

2696850. 

2702087. (0.2) 

2696558. (- 

-0.0) 

28 

2p3p 

^D3 

2107390. 

2110622. ( 

0.2) 

2106614. (-0.0) 

77 

2p4d 

3po 

- 

2702915. ( 

-) 

2697461. ( 

-) 

29 

2p3p 


2119690. 

2123472. ( 

0.2) 

2118706. (-0.0) 

78 

2p 4d 


- 

2705726. ( 

-) 

2700008. ( 

-) 

30 

2p3p 

'"Po 

- 

2132252. ( 

-) 

2126478. ( -) 

79 

2p 4d 


- 

2706671. ( 

-) 

2701074. ( 

-) 

31 

2p3p 

^Pl 

2128680. 

2133785. ( 

0.2) 

2128247. (-0.0) 

80 

2p4f 

‘F3 

- 

2707963. ( 

-) 

2701789. ( 

-) 

32 

2p3p 

^P2 

2130410. 

2135720. ( 

0.2) 

2130292. (-0.0) 

81 

2p4f 

3 F 3 

- 

2708388. ( 

-) 

2702400. ( 

-) 

33 

2p3d 

3 po 

- 

2139901. ( 

-) 

2135677. ( -) 

82 

2p 4d 

3 d° 

- 

2708430. ( 

-) 

2703108. ( 

-) 

34 

2p3d 

3 po 

^3 

- 

2142598. ( 

-) 

2138636. ( -) 

83 

2p4f 

3 F 2 

- 

2708468. ( 

-) 

2702641. ( 

-) 

35 

2p3d 

Ido 

2141580. 

2144437. ( 

0.1) 

2140459. (-0.1) 

84 

2p4f 

^F 4 

- 

2708819. ( 

-) 

2702505. ( 

-) 

36 

2p3d 

3 po 

^ 4 

- 

2145366. ( 

-) 

2141681. ( -) 

85 

2p 4d 

3 po 

- 

2710818. ( 

-) 

2705109. ( 

-) 

37 

2p3p 

‘D2 

2148410. 

2157057.( 

0.4) 

2150745. ( 0.1) 

86 

2p 4d 

3 po 

^1 

- 

2711544. ( 

-) 

2705747. ( 

-) 

38 

2p3d 

^0° 

2160650. 

2165232. ( 

0.2) 

2159945. (-0.0) 

87 

2p 4d 

3 po 

0 

- 

2711943. ( 

-) 

2706094. ( 

-) 

39 

2p3d 

^0° 

2161960. 

2165961. ( 

0.2) 

2160762. (-0.1) 

88 

2p4f 

^P 3 

- 

2712961. ( 

-) 

2706897. ( 

-) 

40 

2p3d 

^D° 

2163340. 

2167469. ( 

0.2) 

2162400. (-0.0) 

89 

2p4f 

^P 4 

- 

2713531. ( 

-) 

2707350. ( 

-) 

41 

2p 3d 

3 po 

2170190. 

2174402. ( 

0.2) 

2169739. (-0.0) 

90 

2p4f 

'P5 

- 

2715225. ( 

-) 

2708632. ( 

-) 

42 

2p3d 

3 po 

^1 

2171680. 

2175473. ( 

0.2) 

2170906. (-0.0) 

91 

2p4f 

3 D 3 

- 

2716473. ( 

-) 

2710743. ( 

-) 

43 

2p3d 

3 po 

0 

- 

2176019. ( 

-) 

2171512. ( -) 

92 

2p4f 

‘D4 

- 

2717105. ( 

-) 

2710218. ( 

-) 

44 

2p3p 

‘So 

- 

2193603. ( 

-) 

2186094. ( -) 

93 

2p4f 

3 D 2 

- 

2717141. ( 

-) 

2711361. ( 

-) 

45 

2p3d 

Ipo 

^3 

2192860. 

2203326. ( 

0.5) 

2195894. ( 0.1) 

94 

2p4p 

‘So 

- 

2717922. ( 

-) 

2706869. ( 

-) 

46 

2p3d 

Ipo 
^ 1 

- 

2208387. ( 

-) 

2201414. ( -) 

95 

2p4f 

^Di 

- 

2718397. ( 

-) 

2712736. ( 

-) 

47 

2s 4s 

“Si 

- 

2477026. ( 

-) 

2472580. ( -) 

96 

2p4f 

‘D2 

- 

2719895. ( 

-) 

2713991.( 

-) 

48 

2s 4s 

‘So 

- 

2488107.( 

-) 

2483578.( -) 

97 

2p 4d 

Ipo 

- 

2722212.( 

-) 

2715080. ( 

-) 

49 

2s 4p 

3 po 

^0 

- 

2504555.( 

-) 

2500383. ( -) 

98 

2p 4d 

Ipo 

^1 

- 

2723448.( 

-) 

2716711.( 

-) 


4.1 Peak abundance temperature 

In Fig. 12 we compare our 98-level CC ICFT and Breit-Pauli R- 
matrix Maxwellian effective collision strengths, T, for all inelastic 
transitions in Al at the temperature of peak abundance, 10'’ K, in 
an electron collisional plasma. We recall that both calculations use 
the exact same (98-level Cl) structure. We note excellent agreement 
between the two scattering methods, as illustrated by the points ly¬ 
ing on the diagonal, only 82 points out of 4753 have a difference 


larger than 10%. This is to be contrasted with the spread of points 
shown in Fig. [T] which results from our two different atomic struc¬ 
tures. 

Only two points differ from the diagonal more than 50%, at 
around T = 10“*, they correspond to the transitions 67 - 75: 
2p4p^D2 - 2p4d *D° and 83 - 98: 2p4f ^F 2 - 2p4d *Pj. Both are 
dipole allowed through spin-orbit mixing only. It is likely that in¬ 
cluding the spin-orbit interaction in the collision calculation causes 
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Figure 2. A comparison of effective collision strengths from 98-level CI/CC 
ICFT versus Breit-Pauli /?-matrix calculations for all inelastic transitions 
amongst the 98 levels of A1 at Te = 10® K. o: transitions with upper level 
n = 2; □: transitions with upper level n = 3; «: transitions with upper level 
« = 4; dashed lines: 20% fractional difference. 



Figure 3. A comparison of effective collision strengths from 98- vs 238- 
level CI/CC ICFT S-matrix calculations for all inelastic transitions amongst 
the 98 lowest common levels of A1 at Te = 10® K. o: transitions with 
upper level n = 2; □: transitions with upper level n = 3; «: transitions with 
upper level n = 4; dashed lines: 20% fractional difference. 

the difference, providing additional mixing via the (A^-l- l)-electron 
Hamiltonian diagonalization. 

Next, in Fig. [3 we compare our 98- versus 238-level CI/CC 
ICFT effective collision strengths as a whole. In contrast to Fig. (2] 
we see a much wider spread of points, i.e. the agreement between 
different scattering methods is much (much) better than that ob¬ 
tained using different configuration interaction basis sets. Like the 
comparison of line strengths and O™® shown in Fig. [T] we note 
the excellent and very good agreement between the two sets of re¬ 
sults for transitions up to n = 2 and n = 3, respectively, while 
there is a much wider spread in the comparison for transitions up 
to n = 4. However, unlike the atomic structure comparison, there 
is a systematic shift above the diagonal. The 238-level CC results 
are systematically larger than the 98-level ones. Thus, we conclude 
that this is not due to the differences in atomic structure, which were 
evenly distributed above and below the diagonal, rather that this is a 
measure of the lack of convergence of the close-coupling expansion 


Table 3. Number of transitions in Figs[T]-|T]which differ by more than a 
certain relative en'or 6 = ITgg — T 238 I/T 238 (or T —* ya,}, as a percentage. 


Rel. eiTor (%) 

Fig.H i/00 

98 vs 238 Cl 

Fig.E] 98 CI/CC T 
BP vs ICFT 

Fig.|3] ICFT T 
98 vs 238 CI/CC 

1 

3778 

1336 

4579 

2 

3600 

803 

4400 

3 

3416 

500 

4243 

4 

3266 

350 

4077 

5 

3141 

260 

3928 

6 

3022 

206 

3798 

7 

2914 

158 

3676 

8 

2804 

127 

3569 

9 

2722 

106 

3460 

10 

2644 

82 

3357 

20 

2068 

22 

2582 

30 

1643 

9 

2090 

40 

1356 

4 

1725 

50 

1163 

2 

1449 

75 

846 

2 

1113 

100 

707 

2 

901 

150 

538 

1 

647 

200 

443 

0 

505 

300 

336 

0 

320 

1000 

187 

0 

88 

Total 

4035 

4753 

4753 


in the 98-level CC calculation for the transitions involving « = 4 
levels. This mirrors the lack of convergence of the 98-level config¬ 
uration interaction expansion for n = 4 levels. We note that there 
are still more than one hundred levels which lie above the n = 4 
levels in the 238-level CI/CC case. 

Specifically, one half of the transitions in Fig.[3]differ by more 
than 30%. The number of transitions differing by more than 20% 
(lines indicated in the plot) correspond to 1 for u = 2 (out of 45), 
256 for n = 3 (out of 990) and 2331 for n = 4 (out of 3718). Note 
also that the 3 transitions which were circled in Fig. [T] are again 
circled in Fig. The strong transitions 75 - 80 and 73-81 do 
illustrate how outliers in the atomic structure comparison (factor 20 
and 30 difference) show-up as outliers in the collision comparison 
(factor 10 and 20 difference). For the weaker transition 31-41 
the resonances in the collision calculation “dampen” the difference 
in atomic structure, it being “just” a factor of 1.6 now. In Table 
we give the exact number of transitions which have a difference 
6 = IT 98 - T 238 I/T 238 larger than a given percentage for the y„ 
for the 98- vs 238-level Cl atomic structures as well as the T at 
r = 10® K for the 98-level CC ICFT vs Breit-Pauli and the 98- vs 
238-level CC ICFT comparisons. 

We can confirm further that the systematic increase of the 
effective collision strengths to u = 4 in the 238-level CC cal¬ 
culation over those of the 98-level CC calculation is due to the 
lack of convergence of the close-coupling expansion in the latter. 
Like the convergence of the Cl expansion, the convergence of the 
close-coupling expansion is essentially independent of the coupling 
scheme, i.e. the specific /^-matrix method used, be it LS, ICFT, 
Breit-Pauli or DARC. We illustrate this in Fig. |4] where we make 
a similar comparison of 54- vs 130-term CC LS -coupling R-matrix 
effective collision strengths. We see the same systematic increase 
for transitions to « = 4 as in the comparison of ICFT R-matrix 
effective collision strengths. 

The plots we have shown so far are useful in the respect that 
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Figure 4. A comparison of effective collision strengths from 54- vs 130- 
term CI/CC LS -coupling R-matrix calculations for all inelastic transitions 
amongst the 54 lowest common terms of Al at Te = 10® K. o: transitions 
with upper term n = 2; □: transitions with upper term n = 3; o: transitions 
with upper term n = A\ dashed lines: 20% fractional difference. 


they allow us to make comparisons as a function of the strength of 
a transition — larger differences are acceptable for weaker transi¬ 
tions. On the other hand, because of the wide range of strengths, 
it is necessary for such plots to be logarithmic. If we plot the ra¬ 
tio of results from two different calculations we can make a linear 
comparison. We do this in Fig. |5] for the same comparison as we 
made in Fig. [3 and with respect to the lower level of the transition 
( |5(a)[ ( or the upper one ( |5(b)[ >. Fig. |5(a)| shows that the transitions 
with the largest scatter are those from 2s 3/ (lower level index 11- 
20) and 2s 4/ (lower level index 47-60) up to n = 4, and where we 
have highlighted by symbol / colour all upper levels with the same 
n-value. It bears a very strong resemblance to th e same comparison 
(figure 2) made bv lAggarwal & KeenanI l l2015ah to compare the 98- 
level CI/CC DARC effective collision strengths with the 238-level 
CI/CC ICFT ones, except that they did not differentiate (highlight) 
the different n-values of the upper levels. In contrast, |5(b)| clari¬ 
fies that the differences in the effective collision strengths become 
increasingly larger as the upper levels excited move closer to the 
last one included in the 98-level CC calculation. This in turn means 
that the effective collision strengths to the uppermost levels of the 
238-level CC calculation of iFemandez-Menchero et ^ ( |2014J) are 
increasingly unconverged with respect to the close-coupling expan¬ 
sion. However, based on the present convergence study to n = 2 and 
n = 3 we can expect that their results for transitions up to « = 4 to 
be well converged, but those to n = 5 less so since only a partial set 
of n = 6,7 levels (to 1 = 2) were included in their close-coupling 
expansion. 

While all of the plots shown so far are visually appealing, es¬ 
pecially Fig. |5] none of them give any indication of the number 
of transitions whose quantities differ by any given amount — we 
cannot tell the density of points close to the diagonal (Figs[T}0 or 
sitting at unity (Fig.|5j. We must be wary of such plots misleading 
us as to the level of agreement, as opposed to disagreement. Only a 
table like Table^gives such an answer. 

Thus, these comparison s demonstrate that the observation by 
lAggarwal & Keen^ j2015al ). that the 238-lev el CC effective colli- 
sions strengths of lFernand^-Menchero et al.l (1201 4al) are system¬ 
atically and increa singly larger with higher exci tations than the 98- 
level CC results of lAggarwal & KeenanI (1201 4ch over a wide range 
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Figure 5. The ratio of effective collision strengths from 98- vs 238-level 
CI/CC ICFT R-matrix calculations vs (a) lower level and (b) upper level 
index, for all inelastic transitions amongst the 98 lowest common levels of 
Al at Tc = 10® K. Positive values indicate Tgg > T 238 and negative 
values T 238 > T 98 . o: transitions with upper level n = 2; □: transitions with 
upper level « = 3; o: transitions with upper level « = 4. 


of temperatures is correct, but it is due to the l ack of convergence of 
the CC expansion of the 98-level CC results of lAggarwal & KeenanI 
( l2014ct) . particularly with respect to the n = 4 levels. 

Finally, there is in fact is good accord between comparable 
calculations, viz. 98-level CI/CC, we make such a comparison of 
ICFT vs DARC in Fig. |6] The increasing difference seen as one 
progresses to higher levels is a reflection of the increasing lack of 
convergence in the atomic structure. While both use the same Cl ex¬ 
pansion, there is no reason for both to give the same unconverged 
result. It is interesting to note that the weakest transitions, mostly 
forbidden one s, show better agreement than some of the stronger 
allowed ones. lAggarwal & KeenanI ( 1201 Sail have already made 
a detailed comparison of 98-level CI/CC DARC and 238-CI/CC 
ICFT effective collisions strengths for transitions from the ground 
state. They highlighted several transitions to « = 4 which were 
discrepant, particularly at high temperatures. We consider them in 
detail in Sec. 14. 31 
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Figure 6. A compaiison of eff ective collision strengths from 98-level CI/CC 
ICFT versus Dirac S-matrix jAggarwal & Keenai]ll2014al calculations for 
all inelastic transitions amongst the 98 levels of A1 at = 10* K. o: 
transitions with upper level n = 2; □: transitions with upper level « = 3; o: 
transitions with upper level « = 4; dashed lines: 20% fractional difference. 


Figure 7. A comparison of effective collision strengths from 98- vs 238- 
level CI/CC ICFT S-matrix calculations for all inelastic transitions amongst 
the 98 lowest common levels of A1 at Te = 2 X 10“* K. o: transitions with 
upper level « = 2; □: transitions with upper level n = 3; «: transitions with 
upper level n = 4; dashed lines: 20% fractional difference. 


4.2 Low temperature 

The effective collision strengths that we have presented so far have 
been relevant to the temperature of peak abundance (of A1 in 
an electron collisional plasma, such as found in the solar atmo¬ 
sphere and magnetic fusion devices. In photoionized plasmas the 
same charge state exists at much lower temperatures. The role of 
resonances becomes more important at low temperatures both with 
respec t to their magnitude, resolut i on and , in particular, their po¬ 
sition. iFemandez-Menchero et al] l l2014ah carried-out an exhaus¬ 
tive analysis of the convergence of the effective collision strengths 
with respect to the energy step used to map-out the resonances. 
They calculated the effective collision strengths T by convolution 
of the ordinary collision strengths with a Maxwellian electron 
energy distribution. Then they reduced the energy step by a factor 
one half and re-calculated T. After repeatedly reducing the energy 
step down to a factor one eighth from the original one, the worst 
case transition from the ground state, the 1 - 78, only changed by 
10% relative to the previous value of T, at the lowest temperature 
calculated of 2 x 1 0** K, and by 1% at 10* K. 

We note that lAggarwal & KeenanI (1201 Sal) incorrectly re¬ 
ported the energvste£usedJn_diejesonance_region_b^Femande^ 
Menche ro et al. l l2014ah . As they stated. lFemandez-Menchero et al.l 
) l2014al) used an energy step in the resonance region that scales as 
~ z in the ion charge, not along the sequence. In particular the 
step length used was 6.94 x 10“*z- Ry for Al 4.89 x 10“*z^ Ry 
for Cl'^+, 4.09 X 10^*2-Ry for K‘*+, 3.30 x 10^*2-Ry for Ti‘®+, 
and 2.02 x 1 0~*2^ Ry for Ge^^^. These steps are c omparable to the 
ones used by lAggarwal & KeenanI ( 120121 I^14all3) . being slightly 
finer forAl_^T_and_shghtl^_coarserfor_Ge^^Fernandez-Menchero 
et al. ( l2014al) used this fine mesh only over 2J = I - 23, which 
corresponds to their exchange calculation. They used their coarse 
energy mesh across the resonance region as well for 27 = 25 - 89. 
Such 7-values can only give rise to high-« resonances, which by 
definition are narrow. A simple calculation with autostructure re¬ 
veals that the stronges t resonances have widths < 10~* Ry d. Thus, 
the results provided bv iFemandez-Menchero et ^ ( l2014ar) are con¬ 
verged with respect to the collision energy step and all significant 
resonances are well resolved. Di fferences in low temperature e ffec- 
tive collision strengths between lAggarwal & KeenM ( l2014ch and 


iFernandez-Menchero et al.l (l2014ah can not be ascribed to the reso¬ 
lution of the resonances. 

The largest source of error in the low temperature effective 
collision strengths arises from the inaccuracy in the positioning of 
resonances which lie just above threshold when the temperature (in 
energy units) starts to become comparable in magnitude with the 
uncertainty in position of these resonances. To a first approxima¬ 
tion, this uncertainty in position is given by the difference between 
the calculated and observed values for the energy level to which 
the resonance is attached. In general, the specific level is not know, 
without detailed resonance analysis. Energy level accuracy can be 
improved theoretically via the use of pseudostates or purely ‘ex¬ 
perimentally’ through the use of observed energies or by a com¬ 
bination of theory and observation using term energy corrections. 
Each has its limitations: the use of pseudostates can lead to pseu¬ 
doresonances at higher energ i es wh ile not all levels may be known 
observationally. IStorev et al.l 120141) discuss the various considera¬ 
tions that need to be made in order to calculate accurate data at very 
low temperatures for planetary nebulae, for example. 

Eig.|7]shows a comparison of 98- and 238-level CC ICET ef¬ 
fective collision strengths for transitions between the lowest com¬ 
mon 98 le vels of Al at a photoion ized plasma temperature of 
2 X 10“* K iKallman & Bauti^l200ll) . The same energy grid was 
used in both cases in the resonance region. Differences are much 
larger than the ones seen at the temperature of peak abundance 
in an electron collisional plasma, even for transitions between the 
low-lying levels (n = 2). The increased differences in the effective 
collision strengths is likely due to the position of the resonances, 
especially where the 238-level CC effective collision strengths are 
smaller than the 98-level ones. In addition, the systematic enhance¬ 
ment of the effective collision strengths to n = 4 levels due to res¬ 
onances attached to n > 4 is increased due to the greater relative 
contribution from resonances; i.e. the lack of convergence of the 
close-coupling expansion for these levels in the 98-level calcula¬ 
tion becomes even more significant. However, at low temperatures, 
most modelling applications involve transitions from the ground 
state, and perhaps a metastable. The factor exp{-AEIkT) arising in 
the excitation rate coefficient means that only the lowest few ex¬ 
cited levels are of interest, i.e. transitions within n = 2. 
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The inaccuracy in the position of the resonances makes ef¬ 
fective colli sion strengths from both Femandez-Menchero et al] 
(l2014al) and lAggarwal & KeenanI ||201211201 4alldl increasingly un- 
reliable at low temperatures. For example, if we shift all resonances 
down in energy by 0.002 Ryd (comparable with the accuracy of 
some energy levels) then the effective collision strength for the 
1-4 : 2s^ 'So - 2s2p^P° transition changes by a factor of 2 at 
2 X 10“* K, although this is reduced to a 20% effect by 2 x 10^ K. 
Nevertheless, such results can be used for estimation purposes. 


4.3 High temperature 


We calculate effective collision strengths over a wide range of tem¬ 
peratures, as defined by the OPEN-ADAS adf04 file format viz. 
2 X 10^ - 2 X lO* (z + 1)^ K, to cover all possible applications 
to electron-collisional plasmas. Ffigher temperature effective col¬ 
lision strengths require ordinary collision strengths to higher en¬ 
ergies, which in turn require the contribution from higher partial 
waves, and this needs to be handled efficiently and accurately by 

Rmiatrixcalculatmns^^ _ 

Aggarwal & Keenan ( 1201 Sail highlighted several transitions 
(1 - 64, 1-70 and 1 - 80) for which they observed large differences 
betweenJhe238devd^O/CCJCFTjesults_of Femandez-Menchero 
et al. i2014ah and the 98devel Cl/CC PARC ones of Aggarwal & 
Keenan~ T2014ch , particularly at high temp eratures. They suggested 
that the use of the lBurgess & Tulhi jl992ll formulae at high energy 
was perhaps a major source of error, i.e. that the ICFT calcula¬ 


tions did not go high enough in ener gy for the collision strengths to 
have reached their asymptotic form. lAggarwal & Keen^ ( 1201 Sail 
also commented-oiMheneglectofelectron^exchangeb^Femandez- 
Menchero et al. j2014al) at high-7. 

The ‘to p-up’ procedure used for angular momentum is 
described in Femandez-Menchero et al.l ( 12014al) . In addition, 
iFemandez-Menchero et al.l ( l2014ah included electron exchange for 
angular momenta up to 27 = 23, and then used a non-exchange 
calculati on for the rest of the angula r momenta calculated: 27 = 
25 - 89. dAggarwal & Keenan ( l2014ch included exchange for all of 
th e angular momenta calculated, up to 27 = 91.) The method used 
bv iFemandez-Menchero et ^ d2014ah is not a source of significant 
inaccuracy. By 27 = 25 the smallest exchange multipole is larger 
than 10. Neglect of higher exchange multipoles causes a small un¬ 
derestimate at the highest temperatures for a few very weak highly 
forbidden transitions, i.e. ones that not only have no target mix¬ 
ing with allowed transitions (i.e. zero limit points) but also are not 
strongly enhanced by coupling. By extending the inclusion of ex¬ 
change to 27 = 51 we find no transition differing by more than 5% 
up to 2 X 10’ K, rising to 10% at 2 x 10* K, for Al 

With regards to energy, IFemandez-Menchero et all ( l2014all 
extended the outer region R-matrix calculation up to three times the 
ionization potential, 88 Ry in the case of Al They then carried- 
out a linear interpolation of the reduced collision strength, y, as a 
function of the reduced scattering energy, x e [ 0 , 1 ], for dipole and 
Bom allowed transitions (while forbidden transitions were extrap¬ 
olated) as follows. 

For a given excitation, let E denote the final scattered energy 
(with EE the excitation energy still) and define 


B = 


E 

AE ■ 


(5) 


Then, at threshold (£ = 0) s = 0. Following iBurgess & Tullvl 
( Il992h . we divide all transitions into one of three cases to represent 




Figure 8. Reduced collision strengths versus energy from 238-level CI/CC 
ICFT R-matrix calculations for selected transitions of Al o: denotes the 
last finite calculated and infinite energy points; dashed line: indicates ex¬ 
trapolation for a forbidden transition. (C = 2.) 


(x,y), based-on their infinite energy values yoo, or lack thereof, as 
described in sec. [3] 


Dipole transitions 
InC 


= 1 


ln(e + C) 


n 


y = 


ln(£ + e) 

Born transitions 
e 


B + C 

= Q. 

• Forbidden transitions 

B 


B + C 
= (e -f 1)“ 


Q. 


( 6 ) 


(7) 


( 8 ) 


where C is an arbitrary visual scaling parameter. In the last case, 
formally, o- = 2 in the infinite energy limit. At high but finite en¬ 
ergies, a more accurate approach is to determine a from two rea¬ 
sonably well-separated high-energy points so as to take account 
of enhancement or depletion of these weak high-energy collision 
strengths by continuum coupling. This is carried-out automatically, 
but restricted to the range a = 1-3, i.e. within reasonable physical 
bounds. 

Fig. shows the reduced collision strength in a Burgess- 
Tully (x, y) plot from the 238-level C I/CC ICFT calculation for 
Al by IFemandez-Menchero et ^ ( l2014all . Fig. | 8 (a)| is for a 
strong dipole transition 1-52 : 2s’ 'So - 2s4p 'Pj, Fig. | 8 (b)| the 
dipole transition 1 - 64 : 2s’'So - 2p4s'Pj which takes place 
through configuration mixing, Fig. | 8 (c)| the Born-allowed transi¬ 
tion 1 -70 : 2 s’ *So - 2p4p ’Po which also takes place only through 
configuration mixing, and Fig. | 8 (d)| is for the forbidden transition 
1 - 80 : 2 s’' So - 2p4f 'F 3 which is a very weak two-electron jump. 
An automatically determined value of a = 1 for this transition was 
used to extrapolate the reduced collision strength as a function of 
reduced energy — see equation]^ Fig.[ 8 ]shows that all of the tran¬ 
sitions have reached the assumed asymptotic form. What then is 
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the source of the d ifferences in high temperature effective collision 
strengths noted bv lAggarwal & KeenanI (l2015all ? The answer lies 
in the atomic structure. 

In Table |4] we compare effective collision strengths from the 
98- and 238- level CI/CC ICFT calculation s with the 98-level CI/CC 
DARC ones jAggarwal & Keenanl2014(J) . We note first that results 
for the strong dipole transition 1-21 are independent of atomic 
structure (98- vs 238-level Cl !/„ = 4S/3) and close-coupling ex¬ 
pansion (98- vs 238-level CC T). However, if we consider the weak 
dipole transition 1 - 64 we see that the limit value is a fac¬ 
tor 6.2 larger for the 98- vs 238-level Cl (Breit-Pauli) case and 
this leads to a factor of 2.08 in the corresponding ICFT effective 
collision strengths at 2 x 10^ K. Indeed, the difference in effective 
collision strengths would likely be larger were it not for the fact 
that the 238-level CI/CC is (much) larger at (much) lower tempera¬ 
tures due to additional reso nances and coupling. The DA RC struc¬ 
ture limit point reported bv lAggarwal & KeenanI ( l2015ah is similar 
to the 98-level Cl Breit-Pauli one. Correspondingly, the 98-level 
CI/CC DARC and ICFT effective collision strengths agree to within 
~ 20% over the entire temperature range shown in Table|4] 

For the case of the (weak) Bom-allowed transition 1 - 70 we 
see a similar trend in the comparisons, viz. differences in structure 
(i/oo) leading to corresponding differences in high temperature ef¬ 
fective collision strengths, strong resonance enhancement at lower 
temperature for the 238- vs 98-level CI/CC results and agreement 
to within 30% between the DARC and ICFT 98-level CI/CC re¬ 
sults. dAggarwal & KeenanI ( l2015all do not report Born limits for 
this transition, but clearly the sensitivity to atomic stmcture we see 
reflected in the two Breit-Pauli Bom limits accounts for the remain¬ 
ing difference.) Finally, for the weak forbidden 1-80 transition 
we note a very similar set of comparisons as for the 1-70 tran¬ 
sition, indeed, the DARC and ICFT 98-level CI/CC results agree 
more closely (20%). The 238-level CI/CC ICFT result is increas¬ 
ingly enhanced by resonances over the 98-level CI/CC one at low 
temperatures, by a factor 15 at 2 x lO"* K. We note that including 
higher-/ exchange multipoles does not change the effective colli¬ 
sion strength to 3 s.f. even at the highest temperature considered. 
There are several equally forbidden transitions (1 - 61,68,69, not 
shown) for which the pattern of dis/agreement is very similar to that 
for the 1 - 80, in all cases. 

I n summary, the results of the 9 8-level CI/CC DARC calcula¬ 
tion of lAggarwal & KeenanI ( 12014ch are in much closer agreement 
(indeed, no significant differences) with the present 98-level CI/CC 
ICFT results than the 238-level CI/CC ones across a wide range 
of temperaturesfor^l_ofjheJransiImnsJii 2 Ui^hted_^Aggarwal 
& Keenan l2015a[) . However, the results of the calculations ob¬ 
tained using the 238-level Cl target have a better converged atomic 
stmcture and, correspondingly, give more accurate effective colli¬ 
sion strengths, especially at high temperatures, while the much bet¬ 
ter convergence of the 238-level close-coupling expansion provides 
more accurate results across a wide range of temperatures. 


5 CONCLUSIONS 

Reliable and accurate electron-impact excitation data are key to the 
successful spectroscopic diagnostic modelling of non-LTE plas¬ 
mas. We have compared and contrasted differences in such data 
for the benchmark Be-like A1 ion which we have calculated us¬ 
ing the R-matrix method. Such differences arise through: 1) differ¬ 
ing approximations of relativistic effects, 2) uncertainties in atomic 
structure and 3) errors due to the lack of convergence of the close¬ 


coupling expansion. Error 3) is quantifiable and can be reduced sys¬ 
tematically and reliably — we illustrated this by comparing new 
98-level and previous 238-level CC ICFT R-matrix calculations. 
We find that effective collision strengths to n = 4 levels are sig¬ 
nificantly enhanced over a wide range of temperatures by coupling 
to n > 4 levels. Uncertainty 2) is quantifiable but is more difficult 
to reduce and constrain as an error — we compared 98-level and 
238-level configuration interaction expansion calculations of line 
strengths and infinite energy plane-wave Bom collision strengths 
to illustrate this point. Again, transitions to n = 4 levels are most 
susceptible to lack of convergence but now of the configuration in¬ 
teraction expansion. Differences 1) between ICFT and Breit-Pauli 
R-matrix treatments of relativistic effects are small, and negligible 
relative to 2) and 3), as is to be expected for an element which 
lies below Zn. We illustrated this by a comparison of new 98-level 
CI/CC ICFT and Breit-Pauli R-matrix effective collision strengths 
which use the exact same atomic structure. We also find good ac¬ 
cord between our 98-level CI/CC results and previous ones from a 
98-level CI/CC Dirac-Coulomb R-matrix calculation, particularly 
for transitions from the ground-level. 

However, based-upon the study of effects 1), 2) and 3), we 
conclude that the original 238-level CI/CC ICFT R-matrix results 
are the most complete to-date with respect to convergence of both 
the configuration interaction and close-coupling expansions and a 
reliable treatment of relativistic effects . Or to put it more simply, 
the earlier 238-level CI/CC ICFT work jFemandez-Menchero et al.l 
l20I4ah has a superiortMgeyo^e98develCI/CC^^RC^onejAg- 
garwal & Keenan T20I4cll and provides more accurate atomic data. 

ThuSjWe^ndJo_be_^lseJhe_^ecent_coniecture_^Aggarwal 
& Keenan ( l20I5ah that the ICFT approach may not be completely 
robust. Their conjecture was based upon a comparison of 98-level 
CI/CC DiracR;matrix_effective_collision_stren 2 th^Aggarwal & 
Keenan lrOIdch with those from 238-ley el CI/CC ICFT R-matrix 
calculationsjFemandez;Menchero_et_ay|20]Aa|^j^nthe^Aggarwal 
& Keenan ^20I5all have failed to appreciate the size of the effect 
which the lack of convergence, in both the close-coupling and con¬ 
figuration interaction expansions, has on transitions to the higher- 
lying states (n = 4 in the case of a 98-level CI/CC expansion for 
A1 ^^). This can only be quantified by extending the expansions. 

There is nothing special about A1 with regards to the conver¬ 
gence of the close-coupling expansion. Except perhaps at the low¬ 
est charge states, as they discussed, the effective collision strengths 
fo r all ions in the Be-like sequence, f rom B+ to Kr’^^ calculated 
bv iFemandez-Menchero et aH (l2014alFI using the 238 level CI/CC 
expansion are the most complete an d reliable to-date. The use of a 
98 level CI/CC ex pansion for A1 dAggar wal & Keenanll20143l . 
Cl'^^, Ge^^_|jlA22arwal_&Keenai|20]Aa|^_and_Ti^^J_(Aggar- 
wal & Keenan lrOlTl) means that these effective collision strengths 
for transitions up to « = 4 are increasingly an underestimate over 
a wide range of temperatures, including the temperature of peak 
abundance. 

'* Fig. 5 of the paper of lFemandez-Menchero et ^ l2014al) compared their 
effective collision strengths for the transition 2s ^ *So — 2s2p^P | of 
with the correspo nding interpola ted results from iKeenanl jl988li . The fit¬ 
ting coefficients of lKeenanl jl988h were taken in numerical form from CHI¬ 
ANTI v7.1 rather than being transcribed from the original paper. However, 
one of the^oefficients_^asjnissin^^^ignj_asjointe^_out^^Aggarwal & 
Keenan T2015^i ~^hir ^iirTe~conectedinTheliexrTeiease of CHIANTI. 
The calculated results of lFemandez-Menchero et a i] <2014j/ for this tran si- 
tion ai'e now only 50% lai'ger than the interpolated ones lKeenanlll988h at 
the temperature of peak abundance (1.8 x 10^ K), for example. 
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Table 4. Comparison of effective collision strengths, T, at selected finite temperatures, and the reduced quantity ya, at infinite temperature. See Tablel^for the 
transition indices i, j. 


DARC ICFT 


i-j 

T(K) 

98-level CI/CC 


98-level CI/CC 



238-level CI/CC 


2 x lO’^ 

10 ® 

2 x 10’ 

2 x 10“* 

10 ® 

2x lO’ 

oo 

2 x 10“* 

10 ® 

2 x 10’ 

OO 

1-52 

1.09 (-2) 

1.01 (-2) 

3.01 (-2) 

1.14 (-2) 

1.02 (-2) 

3.20 (-2) 

2.40 (-2) 

1.53 (-2) 

9.68 (-3) 

3.22 (-2) 

2.64 (-2) 

1-64 

5.49 (-4) 

1.78 (-4) 

1.70 (-4) 

5.80 (-4) 

1.80 (-4) 

1.96 (-4) 

1.01 (-4) 

3.12 (-3) 

2.94 (-4) 

9.40 (-5) 

1.63 (-5) 

1-70 

2.68 (-5) 

1.26 (-5) 

5.29 (-6) 

3.48 (-5) 

1.11 (-5) 

4.24 (-6) 

2.68 (-6) 

6.99 (-5) 

3.10 (-5) 

2.18 (-5) 

2.18 (-5) 

1-80 

4.52 (-5) 

2.51 (-5) 

8.45 (-6) 

4.63 (-5) 

2.47 (-5) 

7.05 (-6) 

- 

7.04 (-4) 

8.66 (-5) 

1.09 (-5) 

- 


Notes. DARC: lAggarwal & KeenanI l2014cl) : 98-level ICFT: present work; 238-level ICFT: lFemandez-Menchero etalJ i2014ah . A (B) denotes A X 10^. 


There is nothing s pecial about Be-like ions. The 136-level CC 
DARC calculations of lAggarwal & KeenanI ( l2014bh for Fe'^^ are 
shown to he a systematic underestimate com pared to the 197-level 
CC ICFT calculations of iLiang et al.l 120lOh (in addition, Liang et 
al. used a much larger Cl expansion 2985 vs 136 levels) — see Del 
Zanna et al (MNRAS to be submitted) for another detailed anal¬ 
ysis similar to the present paper’s. The convergence of the close¬ 
coupling expansion increasingly affects the accuracy of collision 
data to the highest-lying states in all R-matrix calculations. Like¬ 
wise, the accuracy of the atomic structure becomes more uncertain 
for the most highly-excited states of the configuration interaction 
expansion. 

In general, care must be exercised when comparing collision 
data calculated using different atomic structures and / or close¬ 
coupling expansions lest one draws false conclusions. Finally, 
given Figs [T] and [3] and Table [3] we suggest that it is fanciful to 
assign a single accuracy rating of 20%, say, to an entire collision 
data set. 
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APPENDIX A: C^^ 

Just recently, another paper bv lAggarwal & KeenanI ( 1201 5bh has ap¬ 
peared, this time on Be-like C^^, making muc h the same claims that 
have ju st been refuted, quite generally, above. [Aggarwal & KeenanI 
( l2015bl) have extended their CI/CC expansion up to « = 5 this time, 
but this 166-level ex pansion still falls short of the 238-le vel expan¬ 
sion up to nl = 7d of lFernandez-Menchero et ^ ( l2014ah . This may 
reduce the systematic differences (understimates) of their results 
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Validity of the ICFT R-matrix method 


up ton^_^somewha^_com£aredJo^hose_o£Fernandez-Menchero 
et al. (l2014al) but the low-charge state means that the errors and un¬ 
certainties due to the difference between the two atomic structures 
will he much larger than for Al®^. Arguably, a Breit-Pauli /^-matrix 
with pseudo-states calculation is required for to give a defini¬ 
tive representation of the CI/CC expansion for all level-resolved 
transitions. 
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